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Abstract 

Via partial resolution of Abelian orbifolds we present an algorithm for ex- 
tracting a consistent set of gauge theory data for an arbitrary toric variety 
whose singularity a D-brane probes. As illustrative examples, we tabulate the 
matter content and superpotential for a D-brane living on the toric del Pezzo 
surfaces as well as the zeroth Hirzebruch surface. Moreover, we discuss the 
non-uniqueness of the general problem and present examples of vastly different 
theories whose moduli spaces are described by the same toric data. Our meth- 
ods provide new tools for calculating gauge theories which flow to the same 
universality class in the IR. We shall call it "Toric Duality." 



1 Research supported in part by the CTP and the LNS of MIT and the U.S. Department of Energy 
under cooperative research agreement # DE-FC02-94ER40818. A. H. is also supported by an A. P. 
Sloan Foundation Fellowship, a DOE OJI award and by the NSF under grant no. PHY94-07194. 



1 Introduction 



The study of D-branes as probes of geometry and topology of space-time has by 
now been of wide practice (cf. e.g. flU). In particular, the analysis of the mod- 
uli space of gauge theories, their matter content, superpotential and /3-function, as 
world- volume theories of D-branes sitting at geometrical singularities is still a widely 
pursued topic. Since the pioneering work in H, where the moduli and matter con- 
tent of D-branes probing ALE spaces had been extensively investigated, much work 
ensued. The primary focus on (Abelian) orbifold singularities of the type C 2 /2Z„ 
was quickly generalised using McKay's Correspondence, to arbitrary (non-Abelian) 
orbifold sine ularities C 2 /(T C SU(2)), i.e., to arbitrary ALE spaces, m a. 

Several directions followed. With the realisation || ^| that these singularities 
provide various horizons, H [|[ was quickly generalised to a treatment for arbitrary 
finite subgroups T C SU(N), i.e., to generic Gorenstein singularities, by J7J. The 
case of SU(3) was then promptly studied in 0, ||, |1(J using this technique and a 
generalised McKay- type Correspondence was proposed in || |TTJ. Meanwhile, via 
T-duality transformations, certain orbifold singularities can be mapped to type II 
brane-setups in the fashion of [[Tj|]. The relevant gauge theory data on the world 
volume can thereby be conveniently read from configurations of NS-branes, D-brane 
stacks as well as orientifold planes. For C 2 orbifolds, the A and D series have been 
thus treated [12|, [13]], whereas for C 3 orbifolds, the Abelian case of x TLy has been 
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solved by the brane box models fllj, [15|]. First examples of non-Abelian C orbifolds 
have been addressed in Ul6l as well as recent works in 111 



Thus rests the status of orbifold theories. What we note in particular is that once 
we specify the properties of the orbifold in terms of the algebraic properties of the 
finite group, the gauge theory information is easily extracted. Of course, orbifolds 
are a small subclass of algebro-geometric singularities. This is where we move on to 
toric varieties. Inspired by the linear cx-model approach of which provides a rich 
structure of the moduli space, especially in connexion with various geometrical phases 
of the theory, the programme of utilising toric methods to study the behaviour of the 
gauge theory on D-branes which live on and hence resolve certain singularities was 
initiated in Jl9(l . In this light, toric methods provide a powerful tool for studying the 
moduli space of the gauge theory. In treating the F-flatness and D-flatness conditions 
for the SUSY vacuum in conjunction, these methods show how branches of the moduli 
space and hence phases of the theory may be parametrised by the algebraic equations 
of the toric variety. Recent developments in "brane diamonds," as an extension of 
the brane box rules, have been providing great insight to such a wider class of toric 
singularities, especially the generalised conifold, via blown-up versions of the standard 



brane setups [pO|| . Indeed, with toric techniques much information could be extracted 
as we can actually analytically describe patches of the moduli space. 

Now Abelian orbifolds have toric descriptions and the above methodolgy is thus 
immediately applicable thereto. While bearing in mind that though non-Abelian 
orbifolds have no toric descriptions, a single physical D-brane has been placed on 
various general toric singularities. Partial resolutions of C 3 /(2Z 2 x 7L 2 ), such as the 
conifold and the suspended pinched point have been investigated in |21], |22| and brane 



setups giving the field theory contents are constructed by |23], |25], p| . Groundwork 
for the next family, coming from the toric orbifold C 3 /(2Z3 x Z^), such as the del 
Pezzo surfaces and the zeroth Hirzebruch, has been laid in |p6 |. Essentially, given the 
gauge theory data on the D-brane world volume, the procedure of transforming this 
information (F and D terms) into toric data which parametrises the classical moduli 
space is by now well-established. 

One task is therefore immediately apparent to us: how do we proceed in the re- 
verse direction, i.e., when we probe a toric singularity with a D-brane, how do we know 
the gauge theory on its world-volume? We recall that in the case of orbifold theo- 
ries, @ devised a general method to extract the gauge theory data (matter content, 
superpotential etc.) from the geometry data (the characters of the finite group T), 
and vice versa given the geometry, brane-setups for example, conveniently allow us to 
read out the gauge theory data. The same is not true for toric singularities, and the 
second half of the above bi-directional convenience, namely, a general method which 
allows us to treat the inverse problem of extracting gauge theory data from toric data 
is yet pending, or at least not in circulation. 

The reason for this shortcoming is, as we shall see later, that the problem is highly 
non-unique. It is thus the purpose of this writing to address this inverse problem: 
given the geometry data in terms of a toric diagram, how does one read out (at least 
one) gauge theory data in terms of the matter content and superpotential? We here 
present precisely this algorithm which takes the matrices encoding the singularity 
to the matrices encoding a good gauge theory of the D-brane which probes the said 
singularity. 

The structure of the paper is as follows. In Section 2 we review the procedure of 
proceeding from the gauge theory data to the toric data, while establishing nomen- 
clature. In Subsection 3.1, we demonstrate how to extract the matter content and 
F-terms from the charge matrix of the toric singularity. In Subsection 3.2, we exem- 
plify our algorithm with the well-known suspended pinched point before presenting 
in detail in Subsection 3.3, the general algorithm of how to obtain the gauge theory 
information from the toric data by the method of partial resolutions. In Subsection 
3.4, we show how to integrate back to obtain the actual superpotential once the F- 



flatness equations are extracted from the toric data. Section 4 is then devoted to 
the illustration of our algorithm by tabulating the D-terms and F-terms of D-brane 
world volume theory on the toric del Pezzo surfaces and Hirzebruch zero. We finally 
discuss in Section 5, the non-uniqueness of the inverse problem and provide, through 
the studying of two types of ambiguities, ample examples of rather different gauge 
theories flowing to the same toric data. Discussions and future prospects are dealt 
with in Section 6. 



2 The Forward Procedure: Extracting Toric Data 
From Gauge Theories 

We shall here give a brief review of the procedures involved in going from gauge theory 
data on the D-brane to toric data of the singularity, using primarily the notation and 



concepts from [ I9[ . In the course thereof special attention will be paid on how toric 
diagrams, SUSY fields and linear cr-models weave together. 

A stack of n D-brane probes on algebraic singularities gives rise to SUSY gauge 
theories with product gauge groups resulting from the projection of the U (n) theory 
on the original stack by the geometrical structure of the singularity. For orbifolds 
€ k /T, we can use the structure of the finite group T to fabricate product U (rii) gauge 
groups 0, 0, 0]. For toric singularities, since we have only (Abelian) U(l) toroidal 
actions, we are so far restricted to product U(l) gauge groupsQ. In physical terms, 
we have a single D-brane probe. Extensive work has been done in p6[ |l9j to see how 
the geometrical structure of the variety can be thus probed and how the gauge theory 
moduli may be encoded. The subclass of toric singularities, namely Abelian orbifolds, 
has been investigated to great detail H, [|, ^2], |26| and we shall make liberal usage 
of their properties throughout. 

Now let us consider the world-volume theory on the D-brane probe on a toric 
singularity. Such a theory, as it is a SUSY gauge theory, is characterised by its matter 
content and interactions. The former is specified by quiver diagrams which in turn 
give rise to D-term equations; the latter is given by a superpotential, whose partial 
derivatives with respect to the various fields are the so-called F-term equations. 
F and D-flatness subsequently describe the (classical) moduli space of the theory. 
The basic idea is that the D-term equations together with the Fl-parametres, in 
conjunction with the F-term equations, can be concatenated together into a matrix 
which gives the vectors forming the dual cone of the toric variety which the D-branes 
probe. We summarise the algorithm of obtaining the toric data from the gauge theory 
in the following, and to illuminate our abstraction and notation we will use the simple 



Proposals toward generalisations to D-brane stacks have been made 
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Toric Diagram Quiver Diagram 



Figure 1: The toric diagram for the singularity (D 3 / [7L 2 x 7L 2 ) and the quiver diagram 
for the gauge theory living on a D-brane probing it. We have labelled the nodes of 
the toric diagram by columns of G t and those of the quiver, with the gauge groups 
17(1) {A,B,C,D}- 



example of the Abelian orbifold (D 3 /(2Z 2 x 7L 2 ) as given in Figure |I[ 

1. Quivers and D- Terms: 

(a) The bi-fundamental matter content of the gauge theory can be conve- 
niently encoded into a quiver diagram Q, which is simply the (possibly 
directed) graph whose adjacency matrix is precisely the matrix of 
the bi-fundamentals. In the case of an Abelian orbifoldf] prescribed by the 
group T, this diagram is the McKay Quiver (i.e., for the irreps Ri of T, 
is such that R®Ri = ®ja^Rj for some fundamental representation R). We 
denote the set of nodes as Qo '■— {v} and the set of the edges, Q\ := {a}. 
We let the number of nodes be r; for Abelian orbifolds, r = |T| (and for 
generic orbifolds r is the number of conjugacy classes of T). Also, we let 
the number of edges be m; this number depends on the number of su- 
persymmetries which we have. The adjacency matrix (bi-fundamentals) is 

r 

thus r x r and the gauge group is T7 SU(vjj). For our example of 7L 2 x 7L 2 , 

r = 4, indexed as 4 gauge groups U(1)a x U(1)b x U(1)c x U(1)d cor- 
responding to the 4 nodes, while m = 4 x 3 = 12, corresponding to the 

/o 1 1 1 \ 

12 arrows in Figure |l|. The adjacency matrix for the quiver is I { ° I M . 
Though for such simple examples as Abelian orbifolds and conifolds, brane 
3 This is true for all orbifolds but of course only Abelian ones have known toric description. 



setups and 0] specify the values of Wj as well as completely^, there is 
yet no discussion in the literature of obtaining the matter content and 
gauge group for generic toric varieties in a direct and systematic manner 
and a partial purpose of this note is to present a solution thereof. 

(b) From the r x r adjacency matrix, we construct a so-called rxm incidence 
matrix d for Q; this matrix is defined as d v>a := 6 v>head ( a ) - S V)tail ^ for 
v £ Qo and a G Q\. Because each column of d must contain a 1 , a 
—1 and the rest O's by definition, one row of d is always redundant; this 
physically signifies the elimination of an overall trivial U(l) corresponding 
to the COM motion of the branes. Therefore we delete a row of d to define 
the matrix A of dimensions (r — 1) x m; and we could always extract d 
from A by adding a row so as to force each column to sum to zero. This 
matrix A thus contains almost as much information as Oy and once it is 
specified, the gauge group and matter content are also, with the exception 
that precise adjoints (those charged under the same gauge group factor and 
hence correspond to arrows that join a node to itself) are not manifest. 
For our example the 4 x 12 matrix d is as follows and A is the top 3 rows: 
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-1 



The moment maps, arising in the sympletic-quotient language of the toric 
variety, are simply \i :— d ■ |x(a)| 2 where x(a) are the afline coordinates 
of the C r for the torus (C*) r action. Physically, x(a) are of course the 
bi-fundamentals in chiral multiplets (in our example they are -X"ye{ J 4,B,c,£)} 
as labelled above) and the D-term equations for each £7(1) group is ITS 



A = -e a (X)d<«k(a)| 2 -Ci) 

a 

with d the Fl-parametres. In matrix form we have A • |x(a)| 2 = ( and see 
that D-flatness gives precisely the moment map. These ^-parametres will 
encode the resolution of the toric singularity as we shall shortly see. 

2. Monomials and F- Terms: 

(a) From the super-potential W of the SUSY gauge theory, one can write the 
F-Term equation as the system -g^-W = 0. The remarkable fact is that 



l For arbitrary orbifolds, YJ Willi = \T\ where m are the dimensions of the irreps of T; for Abelian 



we could solve the said system of equations and express the m fields in 
terms of r + 2 parametres Vj which can be summarised by a matrix K. 



Xi 



n 



'3 ' 



i = l,2,..,m; j 



1,2,.., 



r + 2 



(2.1) 



This matrix of dimensions m x (r + 2) is the analogue of A in the 
sense that it encodes the F-terms and superpotential as A encodes the 
D-terms and the matter content. In the language of toric geometry K 
defines a conef] M + : a non-negative linear combination of m vectors Ki 
in an integral lattice 7L r+2 . 
For our example, the superpotential is 

W = XacXcdXda — Xac'XcbXba + XqaXabXbc — XcaXadXdc 
-\-XbdXdcXcb ~ XbdXdaXab ~ XdbXbcXcd-, 

giving us 12 F-term equations and with the manifold of solutions parametris- 
able by 4 + 2 new fields, whereby giving us the 12x6 matrix (we here 
show the transpose thereof, thus the horizontal direction corresponds to 
the original fields Xi and the vertical, Vj): 



L 





X-AC 


Xbd 


X-CA 




X-AB 




X-CD 






X-BC 


Xcb 


Xda 


VI 


= Xac 


1 
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1 
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V 2 


= x BD 
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1 





-1 








-1 














V3 


= X B a 

















1 





1 
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V4 


= X C D 














1 





1 








-1 





-1 


V5 


= Xad 








-1 


-1 
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1 








Vg 


= X C B 
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1 




















1 


1 



For example, the third column reads Xca - 
XbdXcb, which the the F-flatness condition 



v 2 v 5 t> 6 , i.e., XadXca = 
aw The details of ob- 



ax DC =o ■ 



taining W and K from each other are discussed in |0|, |26| and Subsection 
3.4. 

(b) We let T be the space of (integral) vectors dual to K, i.e., K ■ T > for 
all entries; this gives an (r + 2) x c matrix for some positive integer c. 
Geometrically, this is the definition of a dual cone N + composed of vectors 
Ti such that K ■ T > 0. The physical meaning for doing so is that K 
may have negative entries which may give rise to unwanted singularities 



5 We should be careful in this definition. Strictly speaking we have a lattice M = 7L r+2 with its 
dual lattice N = 7L r+2 . Now let there be a set of Z5 + -independent vectors {k{\ £ M and a cone is 
defined to be generated by these vectors as a := {J^t a ih | a« 6 H>o}; Our M + should be M n a. 
In much of the literature M+ is taken to be simply := {X)i a j^ I a i 6 ^>o} in which case 
we must make sure that any lattice point contained in M+ but not in must be counted as an 
independent generator and be added to the set of generators {hi}. After including all such points 
we would have = M+. Throughout our analyses, our cone defined by K as well the dual cone 
T will be constituted by such a complete set of generators. 



and hence we define a new set of c fields Pi (a priori we do not know the 
number c and we present the standard algorithm of finding dual cones in 



the Appendix). Thus we reduce ( |2.1| ) further into 



(2.2) 



whereby giving X t — LLj uj 
7L 2 x Z2 example, c = 9 and 



Ylapl^ 3 13 3a with J2jKijTj a > 0. For our 



(- 



T 



J a 



Pi P2 P3 Pi P5 P6 P7 PS P9 \ 



1 
1 



1 

/ 



These new variables p a are the matter fields in Witten's linear u-model. 
How are these fields charged? We have written r + 2 fields v j in terms of 
c fields p a , and hence need c — (r + 2) relations to reduce the independent 
variables. Such a reduction can be done via the introduction of the new 
gauge group £/(l) c ~ < - r+2 ) acting on the p^s so as to give a new set of D- 
terms. The charges of these fields can be written as Qka- The gauge 
invariance condition of Vi under U{\) c ^^ +2 \ by (|2.2|) , demands that the 
(c — r — 2) x c matrix Q is such that J2 a Tj a Qka = 0. This then defines for 
us our charge matrix Q which is the cokernel of T: 



TQ l = {T ja ){Q ka ) t = 0, 



J 



,r+2: a 



c; 



k 



(c-r-2) 



For our example, the charge matrix is (9 — 4 — 2) x 9 and one choice is 
/o 001-11-10 \ 



Q 



ka 



1 1 
1-110-10 







(d) In the linear a-model language, the F-terms and D-terms can be treated 
in the same footing, i.e., as the D-terms (moment map) of the new fields 
p a ; with the crucial difference being that the former must be set exactly 
to zero[] while the latter are to be resolved by arbitrary Fl-parameters. 
Therefore in addition to finding the charge matrix Q for the new fields 
Pa coming from the original F-terms as done above, we must also find 
the corresponding charge matrix Qd for the Pi coming from the original 
D-terms. We can find Qd in two steps. Firstly, we know the charge 



6 Strictly speaking, we could have an F-term set to a non-zero constant. An example of this 
situation could be when there is a term aip + 4>QQ in the superpotential for some chargeless field <fi 
and charged fields Q and Q. The F-term for <f> reads QQ — —a and not 0. However, in our context 
<f) behaves like an integration constant and for our purposes, F-terms are set exactly to zero. 



matrix for Xj under U(l) r 1 , which is A. By ( f2.1| ), we transform the 
charges to that of the v/s, by introducing an (r — 1) x (r + 2) matrix V 
so that V ■ K l 
we have An = 



A. To see this, let the charges of v~ be Vu then by (|2TT| 



V • K t . A convenient V which does so for our 

. Secondly, we use (|2.2| ) 

(4-1) x (4+2) 

to transform the charges from v/s to our final variables p a 's, which is done 



7L 2 x 7L 2 example is 




by introducing an (r + 2) x c matrix Uj a so that U-T l = Id( r+2 ) X (r+2)- 



example, one choice for U is Uj a 



In our 
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oj 



Threfore, combining the two steps, we obtain Qn 

' 1 -1 -1 00 l\ 

-1101000-10 
,-1 01 -1 1 0/ 



(4+2) x 9 

V ■ U and in our 



example, (V ■ U)i a 



3. Thus equipped with the information from the two sides: the F-terms and 
D-terms, and with the two required charge matrices Q and V ■ U obtained, 
finally we concatenate them to give a (c — 3) x c matrix Q t . The trans- 
pose of the kernel of Q t , with (possible repeated columns) gives rise to a 
matrix G t - The columns of this resulting G t then define the vertices of the 
toric diagram describing the polynomial corresponding to the singularity on 
which we initially placed our D-branes. Once again for our example, Qt = 
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-1 

















1 


-1 





-1 
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Ci 
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-1 





C2 


-1 
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-1 
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and G t 





The columns of G t , up to repetition, are precisely marked in the toric diagram 
for 7L 2 x 7L 2 in Figure [l} 

Thus we have gone from the F-terms and the D-terms of the gauge theory to the 
nodes of the toric diagram. In accordance with [27j, Gt gives the algebraic variety 
whose equation is given by the maximal ideal in the polynomial ring 
€[YZ, XYZ, Z, X~ l YZ, XY~ l Z, XZ] (the exponents (i,j,k) in X i Y i Z k are exactly 
the columns), which is uvw = s 2 , upon defining u = (YZ)(XYZ) 2 (Z)(XZ) 2 ;v = 
(YZ) 2 (Z) 2 (X- 1 YZ) 2 ;w = (Zf l (XY '~ l Z)(X Zf and 

s = (YZ) 2 (XYZ)(Z) 2 (X- 1 YZ)(XY- 1 Z)(XZ) 2 ; this is precisely C 3 /(^2 x TL 2 ). In 
physical terms this equation parametrises the moduli space obtained from the F and 
D flatness of the gauge theory 

We remark two issues here. In the case of there being no superpotential we could 
still define i^-matrix. In this case, with there being no F-terms, we simply take K 
to be the identity. This gives T =Id and Q = 0. Furthermore U becomes Id and 



V = A, whereby making Q t = A as expected because all information should now be 
contained in the D-terms. Moreover, we note that the very reason we can construct 
a fT-matrix is that all of the equations in the F-terms we deal with are in the form 
TJ XI 1 = n Xj J ; this holds in general if every field X { appears twice and precisely twice 

i 3 

in the superpotential. More generic situations would so far transcend the limitations 
of toric techniques. 

Schematically our procedure presented above at length, what it means is as fol- 
lows: we begin with two pieces of physical data: (1) matrix d from the quiver encoding 
the gauge groups and D-terms and (2) matrix K encoding the F-term equations. From 
these we extract the matrix G t containing the toric data by the flow-chart: 

Quiver — > d — > A 

I 

F- Terms — * K V *L^ A y 

i i 
T = Dual(fT) U V A U -> 

I 

Q = [Ker(T)]< — Q t = I ^ 1 ^ G t = [Ker(Q t )Y 



3 The Inverse Procedure: Extracting Gauge The- 
ory Information from Toric Data 

As outlined above we see that wherever possible, the gauge theory of a D-brane probe 
on certain singularities such as Abelian orbifolds, conifolds, etc., can be conveniently 
encoded into the matrix Q t which essentially concatenates the information contained 
in the D-terms and F-terms of the original gauge theory. The cokernel of this matrix 
is then a list of vectors which prescribes the toric diagram corresponding to the 
singularity. It is natural to question ourselves whether the converse could be done, 
i.e., whether given an arbitrary singularity which affords a toric description, we could 
obtain the gauge theory living on the D-brane which probes the said singularity. This 
is the inverse problem we projected to solve in the introduction. 




3.1 Quiver Diagrams and F-terms from Toric Diagrams 

Our result must be two- fold: first, we must be able to extract the D-terms, or in 
other words the quiver diagram which then gives the gauge group and matter content; 
second, we must extract the F-terms, which we can subsequently integrate back to 
give the superpotential. These two pieces of data then suffice to specify the gauge 



theory. Essentially we wish to trace the arrows in the above flow-chart from G t back 
to A and K. The general methodology seems straightforward: 



1. Read the column- vectors describing the nodes of the given toric diagram, repeat 
the appropriate columns to obtain G t and then set Q t = Coker(Gi); 

2. Separate the D-term (V ■ U) and F-term (Q t ) portions from Q t ; 

3. From the definition of Q, we obtain]] T = ker(Q). 



4. Farka's Theorem |27| guarantees that the dual of a convex polytope remains 
convex whence we could invert and have K = Dual(T*); Moreover the dual- 
ity theorem gives that Dual(Dual(K)) = K, thereby facilitating the inverse 
procedure. 

5. Definitions U ■ T* = Id and V ■ K* = A (V ■ U) ■ (T* • K*) = A. 

We see therefore that once the appropriate Qt has been found, the relations 

K = Dual(T') A = (V ■ U) ■ (T l ■ K 1 ) (3.3) 

retrieve our desired K and A. The only setback of course is that the appropriate Q t 
is NOT usually found. Two ambiguities are immediately apparent to us: (A) In step 
1 above, there is really no way to know a priori which of the vectors we should repeat 
when writing into the Gt matrix; (B) In step 2, to separate the D-terms and the 
F-terms, i.e., which rows constitute Q and which constitute V ■ U within Q t , seems 
arbitrary. We shall in the last section discuss these ambiguities in more detail and 
actually perceive it to be a matter of interest. Meanwhile, in light thereof, we must 
find an alternative, to find a canonical method which avoids such ambiguities and 
gives us a consistent gauge theory which has such well-behaved properties as having 
only bi-fundamentals etc.; this is where we appeal to partial resolutions. 

Another reason for this canonical method is compelling. The astute reader may 
question as to how could we guarantee, in our mathematical excursion of performing 
the inverse procedure, that the gauge theory we obtain at the end of the day is one 
that still lives on the world-volume of a D-brane probe? Indeed, if we naively traced 
back the arrows in the flow-chart, bearing in mind the said ambiguities, we have no a 
fortiori guarantee that we have a brane theory at all. However, the method via partial 
resolution of Abelian orbifolds (which are themselves toric) does give us assurance. 
When we are careful in tuning the Fl-parametres so as to stay inside cone-partitions 
of the space of these parametres (and avoid flop transitions) we do still have the 

7 As mentioned before we must ensure that such a T be chosen with a complete set of 7L+- 
independent generators; 



resulting theory being physical [26]. Essentially this means that with prudence we 
tune the Fl-parametres in the allowed domains from a parent orbifold theory, thereby 
giving a subsector theory which still lives on the D-brane probe and is well-behaved. 
Such tuning we shall practice in the following. 

The virtues of this appeal to resolutions are thus twofold: not only do we avoid 
ambiguities, we are further endowed with physical theories. Let us thereby present 
this canonical mathod. 



3.2 A Canonical Method: 
Orbifolds 



Partial Resolutions of Abelian 



Our programme is standard p6fl : theories on the Abelian orbifold singularity of the 
form <C k /T for T(k, n) = TL n x 7L n x ...7L n (k — 1 times) are well studied. The complete 
information (and in particular the full Q t matrix) for T(k, n) is well known: k = 2 is 
the elliptic model, k — 3, the Brane Box, etc. In the toric context, k = 2 has been 
analysed in great detail by ||, k = 3, n = 2 in e.g. [p3| , ^5| , 24], k = 3, n = 3 in [26 



Now we know that given any toric diagram of dimension k, we can embed it into 
such a T(k, n)-orbifold for some sufficiently large n; and we choose the smallest such 
n which suffices. This embedding is always possible because the toric diagram for 
the latter is the fc-simplex of length n enclosing lattice points and any toric diagram, 
being a collection of lattice points, can be obtained therefrom via deletions of a subset 
of points. This procedure is known torically as partial resolutions of T(k,n). The 
crux of our algorithm is that the deletions in the toric diagram corresponds to the 
turning-on of the Fl-parametres, and which in turn induces a method to determine a 
Qt matrix for our original singularity from that of T(n, k). 

We shall first turn to an illustrative example of the suspended pinched point 
singularity (SPP) and then move on to discuss generalities. The SPP and conifold as 
resolutions of T(3, 2) = 7L 2 x 7L 2 have been extensively studied in |25| . The SPP, given 



by xy 



zw 



can be obtained from the T(3, 2) orbifold, xyz 



w 



by a single IP 1 



blow-up. This is shown torically in Figure 0. Without further ado let us demonstrate 
our procedure. 

1. Embedding into 7L 2 x TL 2 . Given the toric diagram D of SPP, we recognise that 
it can be embedded minimally into the diagram D' of 7L 2 x 7L 2 . Now information 
on D' is readily at hand ||25|| , as presented in the previous section. Let us re- 
capitulate: 
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Figure 2: The toric diagram showing the resolution of the (D 3 / (7L 2 x Z 2 ) singularity to 
the suspended pinch point (SPP). The numbers i at the nodes refer to the i-th column 
of the matrix G t and physically correspond to the fields pi in the linear cx-model. 



number (i.e., these three- vectors are all co-planar) ensures that D' is Calabi- 
Yau IQ. Incidentally, it would be very helpful for one to catalogue the list 
of Qt matrices of T(3,n) for n = 2,3... which would suffice for all local toric 
singularities of Calabi-Yau threefolds. 

In the above definition of Q' t we have included an extra column (0, 0, 0, 0, C2, C3) 
so as to specify that the first three rows of Q' t are F-terms (and hence exactly 
zero) while the last three rows are D-terms (and hence resolved by Fl-parametres 
0,2.3)- We adhere to the notation in |25| and label the columns (linear a- model 
fields) as pi...pg] this is shown in Figure EL 

2. Determining the Fields to Resolve by Tuning (: We note that if we turn on a 
single Fl-parametre we would arrive at the SPP; this is the resolution of D' 
to D. The subtlety is that one may need to eliminate more than merely 
the 7th column as there is more than one field attributed to each node in 
the toric diagram and eliminating column 7 some other columns correspond- 
ing to the adjacent nodes (namely out of 4,6,8 and 9) may also be elimi- 
nated. We need a judicious choice of ( for a consistent blowup. To do so 
we must solve for fields Pi,..^ and tune the ^-parametres such that at least 
Pi acquires non-zero VEV (and whereby resolved). Recalling that the D- 
term equations are actually linear equations in the modulus-squared of the 
fields, we shall henceforth define Xi := \pi\ 2 and consider linear-systems therein. 




yiilliiiii/ 
which is drawn in Figure [I]. The fact that the last row of Gt has the same 



Therefore we perform Gaussian row-reduction on Q' and solve all fields in 
terms of x 7 to give: x = {x h x 2 , x x + C2 + C3, 2xi - X2+ * r -< 1+ & ; 2 x x - x 2 + ( 2 + 

(- 2xi-X2+X7+Cl+C2+2C3 ™ ^'2+^7-Cl-C2 £2+27+Cl±C2 \ 
^3) 2 > '> 2 ' 2 J' 

The nodes far away from p 7 are clearly unaffected by the resolution, thus the 
fields corresponding thereto continue to have zero VEV. This means we solve the 
above set of solutions x once again, setting £5,1,3,2 = 0, with C 12 ,3 being the vari- 
ables, giving upon back substitution, x = {0, 0, 0, X7 ~ C 2 1 ~ C3 , 0, X7+ b + & , x 7 , X7 ~^ + & , 
Z7+ ^ 1 ~'' 3 }. Now we have an arbitrary choice and we set £3 = and x 7 = £i to 
make p 4 and p§ have zero VEV. This makes P6,7,9 our candidate for fields to be 
resolved and seems perfectly reasonable observing Figure 0. The constraint on 
our choice is that all solutions must be > (since the Xj's are VEV-squared) . 

3. Solving for G t '- We are now clear what the resolution requires of us: in order 
to remove node p 7 from D' to give the SPP, we must also resolve 6, 7 and 9. 
Therefore we immediately obtain G t by directly removing the said columns from 
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G t := coker(Qt) = 

the columns of which give the toric diagram D of the SPP, as shown in Figure 0. 

4. Solving for Q t : Now we must perform linear combination on the rows of Q' t to 
obtain Q t so as to force columns 6, 7 and 9 zero. The following constraints must 
be born in mind. Because G t has 6 columns and 3 rows and is in the null space 
of Qt, which itself must have 9 — 3 columns (having eliminated ^6,7,9), we must 
have 6 — 3 = 3 rows for Q t . Also, the row containing d must be eliminated as 
this is precisely our resolution chosen above (we recall that the Fl-parametres 
are such that £2,3 = and are hence unresolved, while £i > and must be 
removed from the D-terms for SPP). 

We systematically proceed. Let there be variables {ai=i,..,6} so that y : = 
J2i a i T °Wi(Qt) is a row of Q t . Then (a) the 6th, 7th and 9th columns of y 
must be set to and moreover (b) with these columns removed y must be in 
the nullspace spanned by the rows of G t . We note of course that since Q' t was 
in the nullspace of G' t initially, that the operation of row-combinations is closed 
within a nullspace, and that the columns to be set to in Q' t to give Q t are 
precisely those removed in G' t to give G t , condition (a) automatically implies 
(b). This condition (a) translates to the equations {ai + a§ = 0, —at + a 2 — aQ = 
0, — a 2 + 04 = 0} which afford the solution a% = —a^; a 2 = a^ = 0. The fact that 
a 4 = is comforting, because it eliminates the row containing C\. We choose 
ai = 1. Furthermore we must keep row 5 as ( 2 is yet unresolved (thereby setting 



= 1). This already gives two of the 3 anticipated rows of Q t : row 5 and row! 
- row6- The remaining row must corresponds to an F-term since we have ex- 
hausted the D-terms, this we choose to be the only remaining variable: = 1. 
Consequently, we arrive at the matrix 



Qt 



5. Obtaining K and A Matrices: The hard work is now done. We now recognise 
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from Q t that Q 
ker(Q) = 
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K l := Dual(T* 
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to have every entry positive semi-definite. Furthermore we recognise from Q t 



Subsequently we obtain T • K 



, which we do observe indeed 



that V ■ U 
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/ ) , whence we obtain at last, using 
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giving us the quiver diagram (included in Figure |3] for reference), matter content 



and gauge group of a D-brane probe on SPP in agreement with fl25| . We shall 
show in the ensuing sections that the superpotential we extract has similar 
accordance. 



3.3 The General Algorithm for the Inverse Problem 

Having indulged ourselves in this illustrative example of the SPP, we proceed to 
outline the general methodology of obtaining the gauge theory data from the toric 
diagram. 

1. Embedding into € k /{7L n f- x : We are given a toric diagram D describing an 
algebraic variety of complex dimension k (usually we are concerned with local 
Calabi-Yau singularities of k = 2, 3 so that branes living thereon give M = 2, 1 
gauge theories). We immediately observe that D could always be embedded 
into D', the toric diagram of the orbifold C fc / \7L n ) k ~ l for some sufficiently large 
integer n. The matrices Q' t and G' t for D' are standard. Moreover we know that 
the matrix Gt for our original variety D must be a submatrix of G' t . Equipped 
with Q' t and G[ our task is to obtain Q t ; and as an additional check we could 
verify that Q t is indeed in the nullspace of G t . 




Figure 3: The quiver diagram showing the matter content of a D-brane probing the 
SPP singularity. We have not marked in the chargeless field <fi (what in a non-Abelian 
theory would become an adjoint) because thus far the toric techniques do not yet know 
how to handle such adjoints. 



2. Determining the Fields to Resolve by Tuning (: Q' t is a k x a matrix^ (because 
D' and D are dimension k) for some a; G' t , being its nullspace, is thus (a — k) x a. 
D is a partial resolution of D' . In the SPP example above, we performed a single 
resolution by turning on one Fl-parametre, generically however, we could turn 
on as many £'s as the embedding permits. Therefore we let Gt be (a — k) x (a — b) 
for some b which depends on the number of resolutions. Subsequently the Q t 
we need is (k — b) x (a — b). 

Now b is determined directly by examining D' and D; it is precisely the number 
of fields p associated to those nodes in D' we wish to eliminate to arrive at D. 
Exactly which b columns are to be eliminated is determined thus: we perform 
Gaussian row-reduction on Q' t so as to solve the k linear-equations in a variables 
x % := \Pi\ 2 i with F-terms set to and D-terms to Fl-parametres. The a variables 
are then expressed in terms of the Q's and the set B of x^s corresponding to 
the nodes which we definitely know will disappear as we resolve D' —>■ D. The 
subtlety is that in eliminating B, some other fields may also acquire non zero 
VEV and be eliminated; mathematically this means that Order(5) < b. 
Now we make a judicious choice of which fields will remain and set them to 
zero and impose this further on the solution Xj = i r . )0 = Xi(Q; B) from above 
until Order (B) = b, i.e., until we have found all the fields we need to eliminate. 
We know this occurs and that our choice was correct when all Xj > with 
those equaling corresponding to fields we do not wish to eliminate as can be 
observed from the toric diagram. If not, we modify our initial choice and repeat 
until satisfaction. This procedure then determines the b columns which we wish 
8 We henceforth understand that there is an extra column of zeroes and £'s. 



to eliminate from Q' t . 



3. Solving for Gt and Qt'- Knowing the fields to eliminate, we must thus perform 
linear combinations on the k rows of Q' t to obtain the k — b rows of Q t based upon 
the two constraints that (1) the b columns must be all reduced to zero (and thus 
the nodes can be removed) and that (2) the k — b rows (with b columns removed) 
are in the nullspace of G t - As mentioned in our SPP example, condition (1) 
guarantees (2) automatically. 

In other words, we need to solve for k variables {£j=i,..,/c} such that 

k 

Y, x i (Q'th = for j = Pi>P2, -Pb e B. (3.4) 

i=l 

Moreover, we immediately obtain Gt by eliminating the b columns from G' t . 

k 

implies that E E %i (Qt)ij (Gt)mj = for 

i=lj¥v\...b 

m — 1, a — k and hence guarantees that the Qt we obtain is in the nullspace 
of G t . 

We could phrase equation ( |3.4j ) for Xi in matrix notation and directly evaluate 

Qt = NullSpace(Wy ■ Q' t (3.5) 

where Q' t is Q' t with the appropriate columns (pi...b) removed and W is the 
matrix constructed from the deleted columns. 

4. Obtaining the K Matrix (F-term): Having obtained the (k — b)x (a — b) matrix 
Qt for the original variety D, we proceed with ease. Reading from the extraneous 
column of Fl-parametres, we recognise matrices Q (corresponding to the rows 
that have zero in the extraneous column) and V ■ U (corresponding to those 
with combinations of the unresolved £'s in the last column). We let V ■ U be 
c x (a — b) whereby making Q of dimension (k — b — c) x (a — b). The number 
c is easily read from the embedding of D into D' as the number of unresolved 
Fl-parametres. 

From Q, we compute the kernel T, a matrix of dimensions (a — b) — (k — b — 
c) x (a — b) = (a — k + c) x (a — b) as well as the matrix K l of dimensions 
(a — k + c) x d describing the dual cone to that spanned by the columns of T. 
The integer d is uniquely determined from the dimensions of T in accordance 
with the algorithm of finding dual cones presented in the Appendix. From these 
two matrices we compute T* • K l , of dimension (a — b) x d. 

5. Obtaining the A Matrix (D-term): Finally, we use (|3~3| ) to compute (V ■ U) ■ 
(T* • K l ), arriving at our desired matrix A of dimensions c x d, the incidence 



Indeed, as discussed earlier, (3A 



matrix of our quiver diagram. The number of gauge groups we have is therefore 
c + 1 and the number of bi-fundamentals, d. 

Of course one may dispute that finding the kernel T of Q is highly non-unique 
as any basis change in the null-space would give an equally valid T. This is 
indeed so. However we note that it is really the combination T* ■ K l that we 
need. This is a dot-product in disguise, and by the very definition of the dual 
cone, this combination remains invariant under basis changes. Therefore this 
step of obtaining the quiver A from the charge matrix Q t is a unique procedure. 

3.4 Obtaining the Superpotential 

Having noticed that the matter content can be conveniently obtained, we proceed to 
address the interactions, i.e., the F-terms, which require a little more care. The matrix 
K which our algorithm extracts encodes the F-term equations and must at least be 
such that they could be integrated back to a single function: the superpotential. 

Reading the possible F-flatness equations from K is ipso facto straight-forward. 
The subtlety exists in how to find the right candidate among many different linear 
relations. As mentioned earlier, K has dimensions m x (r — 2) with m corresponding 
to the fields that will finally manifest in the superpotential, r — 2, the fields that solve 
them according to ( |2.1| ) and (|2.2j ); of course, m > r — 2. Therefore we have r — 2 
vectors in 2Z m , giving generically m — r + 2 linear relations among them. Say we 
have rowi + row3 — row7 = 0, then we simply write down XiX 3 = X 7 as one of the 
candidate F-terms. In general, a relation diKy = with G 7L implies an F-term 

i 

Y[X?* = 1 in accordance with ( |2.1| ). Of course, to find all the linear relations, we 

i 

simply find the Zj-nullspace of K l of dimension m — r + 2. 

Here a great ambiguity exists, as in our previous calculations of nullspaces: any 
linear combinations therewithin may suffice to give a new relation as a candidate F- 
termf]. Thus educated guesses are called for in order to find the set of linear relations 
which may be most conveniently integrated back into the superpotential. Ideally, we 
wish this back-integration procedure to involve no extraneous fields (i.e., integration 
constants^]) other than the m fields which appear in the K-matrix. Indeed, as we 
shall see, this wish may not always be granted and sometimes we must include new 
fields. In this case, the whole moduli space of the gauge theory will be larger than 
the one encoded by our toric data and the new fields parametrise new branches of 

9 Indeed each linear relation gives a possible candidate and we seek the correct ones. For the 
sake of clarity we shall call candidates "relations" and reserve the term "F-term" for a successful 
candidate. 

10 By constants we really mean functions since we are dealing with systems of partial differential 
equations. 



the moduli in the theory. 

Let us return to the SPP example to enlighten ourselves before generalising. We 



recall from subsection 3.2, that K 





X! 


x 2 


x 3 


x 4 


x 5 


x 6 


VI 


1 

















V2 








1 





1 





V3 





1 














V4 








1 


1 








V5 











1 





1 



from which we 



can read out only one relation X 3 X 6 — X 4 X 5 = using the rule described in the 
paragraph above. Of course there can be only one relation because the nullspace of 
K l is of dimension 6 — 5 = 1. 

Next we must calculate the charge under the gauge groups which this term carries. 
We must ensure that the superpotential, being a term in a Lagrangian, be a gauge in- 
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variant, i.e., carries no overall charge under A. From d - 

we find the charge of X 3 X 6 to be (g^, qB, qc) = (0 + 0, 
of course by our very construction, X4X5 has the same charge. Now we have two 
choices: (a) to try to write the superpotential using only the six fields; or (b) to 
include some new field <fi which also has charge (0, 0, 0). For (a) we can try the ansatz 
W = XiX 2 (X 3 X 6 — X4X5) which does give our F-term upon partial derivative with 
respect to X\ or X 2 . However, we would also have a new F-term XiX 2 X 3 = by 
t^t-, which is inconsistent with our K since columns 1, 2 and 3 certainly do not add 
to 0. 

This leaves us with option (b), i.e., W = 4>(X 3 Xq — X 4 A 5 ) say. In this case, when 
cf) = we not only obtain our F-term, we need not even correct the matter content A. 
This branch of the moduli space is that of our original theory. However, when 7^ 0, 
we must have X 3 = X 4 = X 5 = X 6 = 0. Now the D-terms read \Xi\ 2 — \X 2 \ 2 = — Ci = 
( 2 , so the moduli space is: {0 G (C, X\ G C} such that Ci + C2 = for otherwise there 
would be no moduli at all. We see that we obtain another branch of moduli space. As 
remarked before, this is a general phenomenon when we include new fields: the whole 
moduli space will be larger than the one encoded by the toric data. As a check, we 
see that our example is exactly that given in p5| , after the identification with their 
notation, F 12 — * A 6 , X 2A — > X 3 , Z 23 — > X\, Z 32 — > X 2 , Y u — > A 4 , A 13 — > X 5 , Z 41 — > <f> 
and (XiX 2 — 0) — > (p. We note that if we were studying a non-Abelian extension to 
the toric theory, as by brane setups (e.g. p5| ) or by stacks of probes (in progress 
from |26|]), the chargeless field would manifest as an adjoint field thereby modifying 
our quiver diagram. Of course since the study of toric methods in physics is so far 
restricted to product U{1) gauge groups, such complexities do not arise. To avoid 
confusion we shall henceforth mark only the bi-fundamentals in our quiver diagrams 
but will write the chargeless fields explicit in the superpotential. 

Our agreement with the results of is very reassuring. It gives an excellent 



example demonstrating that our canonical resolution technique and the inverse algo- 
rithm do indeed, in response to what was posited earlier, give a theory living on a 
D-brane probing the SPP (T-dual to the setup in However, there is a subtle 

point we would like to mention. There exists an ambiguity in writing the superpoten- 
tial when the chargeless field 4> is involved. Our algorithm gives W = 4>(X^Xq— X^Xq) 
while P3 gives W = (XiX 2 — 0)(X 3 X 6 — X 4 X 6 ). Even though they have identical 



moduli, it is the latter which is used for the brane setup. Indeed, the toric methods 
by definition (in defining A from Oy) do not handle chargeless fields and hence we 
have ambiguities. Fortunately our later examples will not involve such fields. 

The above example of the SPP was a naive one as we need only to accommodate a 
single F-term. We move on to a more complicated example. Suppose we are now given 

/Xi Xi X3 X4 X5 Xq Xj Xg Xcj 
101000100 



Xi X2 X3 X4 X5 Xq X7 Xg Xg Xiq 

. .. -1 -1 1 1 . 

a = \ b 1 -1 0-10 1 anaK 

CO 1 1 1-1-1-1 

DO 1-11-11-10 



11000010 
100100001 
010101000 
,001010100 
\ 000 1100 

We shall see in the next section, that these arise for the del Pezzo 1 surface. Now the 



nullspace of K has dimension 10 — 6 = 4, we could obtain a host of relations from vari- 
ous linear combinations in this space. One relation is obvious: X2X7 — X^X^ = 0. The 
charge it carries is (qA, qB, qc, Qd) = (0 + 0, —1 + 0, + 1, l + (— 1)) = (0, —1, 1, 0) which 
cancels that of Xg. Hence X 9 (X 2 X 7 — X 3 X 6 ) could be a term in W. Now t^- thereof 
gives X 7 Xg and from K we see that X7X9 — X1X5X10 = 0, therefore, — XiX 2 X 5 Xi 
could be another term in W . We repeat this procedure, generating new terms as we 
proceed and introducing new fields where necessary. We are fortunate that in this case 
we can actually reproduce all F-terms without recourse to artificial insertions of new 
fields: W = X2X7X9 — X^X^Xg — X4X8X7 — X1X2X5X10 + X3X4X10 + XiXsXgXg. 

Enlightened by these examples, let us return to some remarks upon generalities. 
Making all the exponents of the fields positive, the F-terms can then be written as 

n*r=n^> ( 3 - 6 ) 

i j 

with di, bj G 7L + . Indeed if we were to have another field X& such that k $ {i}, {j} 

then the term X^ I n^ ai — 11 -XV ) , on the condition that X^ appears only this once, 
V < i ) 

must be an additive term in the superpotential W. This is because the F-flatness 
condition = implies ( |3.6| ) immediately. Of course judicious observations are 
called for to (A) find appropriate relations (|3.6|) and (B) find X^ among our m fields. 
Indeed (B) may not even be possible and new fields may be forced to be introduced, 
whereby making the moduli space of the gauge theory larger than that encodable by 
the toric data. 

In addition, we must ensure that each term in W be chargeless under the product 



gauge groups. What this means for us is that for each of the terms X k I n Xf 1 - U x\ j 



3 
3 



we must have Charge s (X fe ) + J2 ^Charge,, (A;) = for s = 1, ..,r indexing through 

i 

our r gauge group factors (we note that by our very construction, for each gauge 



group, the charges for JJX" 4 an d for n A 6j are equal). If X k in fact cannot be found 



* 3 

among our m fields, it must be introduced as a new field <fi with appropriate charge. 
Therefore with each such relation ( [3.6|) read from K, we iteratively perform this said 
procedure, checking A sk +Y, ^i^si = at each step, until a satisfactory superpotential 

i 

is reached. The right choices throughout demands constant vigilance and astuteness. 

4 An Illustrative Example: the Toric del Pezzo 
Surfaces 

As the C 3 / (JL 2 x 7L 2 ) resolutions were studied in great detail in p5| , we shall use the 
data from to demonstrate the algorithm of finding the gauge theory from toric 
diagrams extensively presented in the previous section. 

The toric diagram of the dual cone of the (parent) quotient singularity € 3 /(7L 3 x 
7Z 3 ) as well as those of its resolution to the three toric del Pezzo surface are presented 
in Figure |j. 

del Pezzo 1: Let us commence our analysis with the first toric del Pezzo surface^. 
From its toric diagram, we see that the minimal 7L n x 7L n toric diagram into which 
it embeds is n = 3. As a reference, the toric diagram for C 3 /(Z 3 x 2Z 3 ) is given m 
Figure ^ and the quiver diagram, given later in the convenient brane-box form, in 
Figure |5|. Luckily, the matrices Q' t and G' t for this Abelian quotient is given in [^6 



n Now some may identify the toric diagram of del Pezzo 1 as given by nodes (using the no- 
tation in Figure |) (1,-1,1), (2,-1,0), (-1,1,1), (0,0,1) and (-1,0,2) instead of the one we 
have chosen in the convention of pfl, with nodes (0,-1,2), (0,0,1), (—1,1,1), (1,0,0) and 
(0,1,0). But of course these two Gt matrices describe the same algebraic variety. The for- 
mer corresponds to Spec (€{XY~ 1 Z, X 2 Y~ X , X~ l YZ, Z, X~ 1 Z 2 ]) while the latter corresponds 
to Spec (<D[Y- 1 Z 2 ,Z,X- 1 YZ,X,Y}). The observation that (X 2 Y~ r ) = (X)(X- 1 YZ)- 1 {Z), 
(XY- X Z) = (X)(Y)- 1 (Z) and {X^Z 2 ) = {Y' 1 Z 2 )^)^- 1 ) for the generators of the polyno- 
mial ring gives the equivalence. In other words, there is an SL(5, 7L) transformation between the 5 
nodes of the two toric diagrams. 



Adding the extra column of Fl-parametres we present these matrices belowQ: 
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12 In | p6| , a canonical ordering was used; for our purposes we need not belabour this point and use 
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According to our algorithm, we must perform Gaussian row-reduction on Q' t to 
solve for 42 variables Xj. When this is done we find that we can in fact express all 
variables in terms of 3 XiS together with the 8 Fl-parametres Q. We choose these 
three x^s to be £10,29,36 corresponding to the 3 outer vertices which we know must be 



37 



38 



del Pezzo 



Resolving 
10, 36, 29, 38, 32 
5,20,9, 11,26,4, 16, 23 
27, 28, 41, 42, 39, 40 
1,2,3,6, 12, 15, 18, 19 
21,22,24,25,30,31,33,34,35 



Resolving 
10, 13,20, 11,26, 16, 23 
27,28, 36,41,42, 39,40,29 
1,2,3, 17, 19,21,22,24,25 
31,32, 33, 34, 35 



Resolving 
10,5, 13,20,21, 11,26 
16, 23, 27, 28, 41, 42 
38, 39,40,29, 1,2,3 
6, 12, 17, 19,22,24,25 
20,31,32,33,34,3536 




del Pezzo 3 




Resolving 
10, 5, 13, 20, 11,26 
16, 23,27,28,41,42 
38, 39,40,29, 1,2,3 

12, 17, 19,22,24,25 
20,31,32, 33, 34,35 



Hirzebruch 
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Resolving 
10, 36, 29 

5,20, 9, 11, 26,4, 16, 23 
27, 28, 41, 42, 39, 40 
1,2, 3, 6, 12, 15, 18, 19 
21,22,24,25,30,31,33,34,35 
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Figure 4: The resolution of the Gorenstein sine; ularity C 3 /(^3 x 2Z 3 ) to the three 
toric del Pezzo surfaces as well as the zeroth Hirzebruch surface. We have labelled 
explicitly which columns (linear cr-model fields) are to be associated to each node in 
the toric diagrams and especially which columns are to be eliminated (fields acquiring 
non-zero VEV) in the various resolutions. Also, we have labelled the nodes of the 
parent toric diagram with the coordinates as given in the matrix G t for (D 3 / (2Z 3 x 2Z 3 ). 



resolved in going from C 3 /(^3 x 2Z 3 ) to del Pezzo 1. 

Next we select the fields which must be kept and set them to zero in order to 
determine the range for Q. Bearing in mind the toric diagrams from Figure |J these 
fields we judiciously select to be: ^13,8,37,38- Setting £13,8,37,38 = gives us the solution 
{Ce = 0;£ 29 = C7 = C3 = Ci - faxio = d + C5 + (3^36 = C7 - Cs}, which upon 
back-substitution to the solutions £; we obtained from Q' t , gives zero for £13,8,37,33 
(which we have chosen by construction) as well as £7,14,17,32; for all others we obtain 
positive values. This means precisely that all the other fields are to be eliminated 
and these 8 columns { 13, 8, 37, 38, 7,14,17,32 } are to be kept while the remaining 
42-8=34 are to be eliminated from Q' t upon row-reduction to give Qt- In other words, 
we have found our set B to be {1,2,3,4,5,6,9,10,11,12,15,16,18,19,20,21,22,23,24,25, 
26,27,28,29,30,31,33,34,35,36,39,40,41,42} and thus according to (fTBJ) we immediately 
obtain 
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We note of course that 5 out of the 8 Fl-parametres have been eliminated automati- 
cally; this is to be expected since in resolving C 3 / (2L 3 x 2Z 3 ) to del Pezzo 1, we remove 
precisely 5 nodes. Obtaining the D-terms and F-terms is now straight-forward. Using 
( P-3| ) and re-inserting the last row we obtain the D-term equations (incidence matrix) 
to be 
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From this matrix we immediately observe that there are 4 gauge groups, i.e., £7(1) 4 
with 10 matter fields Xi which we have labelled in the matrix above. In an equivalent 
notation we rewrite d as the adjacency matrix of the quiver diagram (see Figure ||) 
for the gauge theory: 




The K-matrix we obtain to be: 



Xi X2 X 3 X4 X 5 Xa X~i Xg Xg Xio • 
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i / 



which indicates that the original 10 fields Xj can be expressed in terms of 6. This was 
actually addressed in the previous section and we rewrite that pleasant superpotential 



here: 



W — X2X7X9 — X^XqXq — ^4X3X7 — X1X2X5X10 + X3X4X10 + XiX^XqXs- 



del Pezzo 2: Having obtained the gauge theory for del Pezzo 1, we now repeat 
the above analysis for del Pezzo 2. Now we have the Fl-parametres restricted as 
{^36 = C2 = 0; C3 = C4> x 29 = C4 + Ce; ^10 = Ci + C4}) making the set to be eliminated 
as B = { 1, 2, 3, 5, 6, 10, 11, 13, 16, 17, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 
32, 33, 34, 35, 38, 39, 40, 41, 42 }. Whence, we obtain 



Qt 



and observe that 4 D-terms have been resolved, as 4 nodes have been eliminated from 
C 3 /(Z^3 x ^3)- From this we easily extract (see Figure ||) 
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moreover, we integrate the F-term matrices 



/ Xl X2 Xz X4 X 5 Xq X-j Xg Xg XlQ Xu X12 X13 \ 



K 



to obtain the superpotential 

W = X 2 XgXu — XgX 3 X 10 — X4X 8 Xn — X 1 X 2 X 7 X 13 + X 13 X 3 X 6 
— X5X12X6 + X1X5X8X10 + X4X7X12. 
del Pezzo 3: Finally, we shall proceed to treat del Pezzo 3. Here we have the range 
of the Fl-parametres to be {Ci = Ce = Ce = 0; x 29 = C3 = -(5] xiq = C4; C2 = ^3e! Cs = 
-C2 - Go}, which gives the set B as {1, 2, 3, 10, 11, 13, 16, 17, 19, 20, 21, 22, 23, 24, 
25, 26, 27, 28, 29, 31, 32, 33, 34, 35, 36, 39, 40, 41, 42}, and thus according to O 
we immediately obtain 
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We note indeed that 3 out of the 8 Fl-parametres have been automatically resolved, 
as we have removed 3 nodes from the toric diagram for € 3 /(ZZ 3 x 2Z 3 ). The matter 
content (see Figure |[) is encoded in 
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and from the F-terms 
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we integrate to obtain the superpotential 

W = XzX$Xi3 — XgXgXn — X5X6X13 — X1X3X4X10X12 
-\-XtXqXi2 + XiX2X^X\qXh + X^XqXi^ — X2X-jX\^. 

Note that we have a quintic term in W; this is an interesting interaction indeed. 

del Pezzo 0: Before proceeding further, let us attempt one more example, viz., the 
degenerate case of the del Pezzo as shown in Figure |j. This time we note that the 
ranges for the Fl-parametres are {C5 = —%29 + C6 — A\ £ 6 = x 29 — B; x 2 g = B + C; ( 8 = 
— £36 + B; 236 = B + C + D; Xio = A + E} for some positive A, B, C, D and E, that 
B = { 1, 2, 3, 4, 5, 6, 9, 10, 11, 12, 15, 16, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 
29, 30, 31, 32, 33, 34, 35, 36, 38, 39, 40, 41, 42 } and whence the charge matrix is 



Qt 
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1 -1 C2+C6+C8 

-1 1 Cl+C3+C5 

-111-1-11 



We extract the matter content (see Figure ||) as d 



and the F-terms as K l 



010 
1-11 

X\ X2 X3 X4 X5 X(i X7 Xs Xg \ 



Xi X2 X3 X4 Xrj X§ X7 X g 

-10-10-10 1 1 1 
1 1 1-1-1-1 

-11-10 



1 





and from the latter 



we integrate to obtain the superpotential 

W — X\X^Xg — X4X5X7 — X2X3X9 



XiXqXq + X2X5X8 + X^XqXy. 



Of course we immediately recognise the matter content (which gives a triangular 
quiver which we shall summarise below in Figure |^) as well as the superpotential 



from equations (4.7-4.14) of ||19|| ; it is simply the theory on the Abelian orbifold 
C 3 /2& 3 with action (a G 2Z 3 ) : {z 1} z 2} z 3 ) — > (e^ Zi, e^f- z 2 , z 3 ). Is our del Pezzo 
then (D 3 /5Z 3 ? We could easily check from the G t matrix (which we recall is obtained 
from G' t of (C 3 /(2Z3 x 7L^) by eliminating the columns corresponding to the set B): 



I -1 1 

G t = [ o i -looo 

V 1 1 2 110 



These columns (up to repeat) correspond to monomials Z,X~ 1 YZ,Y~ 1 Z 2 ,X in the 
polynomial ring <C[X, Y, Z]. Therefore we need to find the spectrum (set of max- 
imal ideals) of the ring (C[Z, X _1 YZ, Y~ X Z 2 , X], which is given by the minimal 
polynomial relation: (X~ l Y Z) ■ (Y~ l Z 2 ) ■ X = (Z) s . This means, upon defining 
p = X~ l YZ\ q = Y~ X Z 2 \ r = X and s — Z, our del Pezzo is described by pqr = s 3 
as an algebraic variety in C 4 (p, q,r, s), which is precisely C 3 /S 3 . Therefore we have 
actually come through a full circle in resolving C 3 /(!Zi 3 x 2Z 3 ) to (D 3 /!Zi 3 and the va- 
lidity of our algorithm survives this consistency check beautifully. Moreover, since 
we know that our gauge theory is exactly the one which lives on a D-brane probe 
on C 3 /Z 3 , this gives a good check for physicality: that our careful tuning of FI- 
parametres via canonical partial resolutions does give a physical D-brane theory at 
the end. We tabulate the matter content and the superpotential W for the del 
Pezzo surfaces below, and the quiver diagrams, in Figure [|. 
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Superpotential W 



X2X7X9 — X-^X^Xo, 
— X^X^Xi — X1X2X5X10 
+X3X4X10 + X\X^X(,X% 



X2X9-X11 — X9X3X10 
— X^X^Xw — X\XiXiX\z 

+X13X3X6 — X5X12X6 
+X1X5X8X10 + X4X7-X12 



XzX%X\z — XgXgXn 

— X^X(,X\z — X1X3X4X10X12 
+^7X9X12 + X1X2X5X10X11 
+X4X6X14 — X2X7X14 





del Pezzo S C 6 /7L a 


Hirzebruch = P 1 X IP 1 


:= F = E x 


Matter a;j 


/0 3 0\ 

3 
\3 ) 






/0 2 2\ 
2 0' 
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\0 2 0/ 




Superpotential W 


X1X4X9 — X4X5X7 
— X2X3X9 — X\X§X% 
+X2X5X8 + X^X^Xi 


XiXgXio — X3X7X10 
— XiX%X§ — X\X§X\2 
+X3X6X11 + X4X7X9 
+X2X5X12 — X4X5JY11 



Upon comparing Figure f| and Figure [5], we notice that as we go from del Pezzo to 
3, the number of points in the toric diagram increases from 4 to 7, and the number of 
gauge groups (nodes in the quiver) increases from 3 to 6. This is consistent with the 
observation for M = 1 theories that the number of gauge groups equals the number 
of perimetre points (e.g., for del Pezzo 1, the four nodes 13, 8, 37 and 38) in the toric 
diagram. Moreover, as discussed in [pB ], the rank of the global symmetry group {Ei 




Hirzebruch 

Figure 5: The quiver diagrams for the matter content of the brane world- volume 
gauge theory on the 4 toric del Pezzo singularities as well as the zeroth Hirzebruch 
surface. We have specifically labelled the U(l) gauge groups (A, B, ..) and the bi- 
fundamentals (1, 2, ..) in accordance with our conventions in presenting the various 
matrices Qt, A and K. As a reference we have also included the quiver for the parent 
22 3 x 2Z 3 theory 



for del Pezzo i) which must exist for these theories equals the number of perimetre 
point minus 3; it would be an intereting check indeed to see how such a symmetry 
manifests itself in the quivers and superpotentials. 

Hirzebruch 0: Let us indulge ourselves with one more example, namely the 0th 
Hirzebruch surface, or simply IP 1 x IP 1 := Fq := E\. The toric diagram is drawn in 
Figure f|. Now the Fl-parametres are {£ 4 = — x 2 g — x 36 — (5 — ($ — A; £ 5 = —A — B\ ( 7 = 
X10 + X29 + %36 + Cs - C; Cs = -X10 - x 29 - x 36 + D;D = A + B;C = A + B;A = 
x\o — E; xiq = E + F; x 2 9 = B + G} for positive A, B, C, D, E, F and G. Moreover, 
B = { 1, 2, 3, 5, 6, 10, 11, 13, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 
32, 33, 34, 35, 36, 38, 39, 40, 41, 42 }. We note that this can be obtained directly by 
partial resolution of fields 21 and 36 from del Pezzo 2 as is consistent with Figure f|. 
Therefrom we obtain the charge matrix 
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from which we have the matter content d = \ -1 -1 1 

01010 
10100 
the quiver for which is presented in Figure [|. The F-terms are 
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from which we obtain 

W = X\X%Xiq— X3X7X10— X2XgXg— X1X6X12+X3X6X11+X4X7X9+X2X5X12— X4X5X11 



a perfectly acceptable superpotential with only cubic interactions. We include these 
results with our table above. 



5 Uniqueness? 

In our foregoing discussion we have constructed in detail an algorithm which calculates 
the matter content encoded by A and superpotential encoded in K, given the toric 
diagram of the singularity which the D-branes probe. As abovementioned, though this 
algorithm gives one solution for the quiver and the i^-matrix once the matrix Q t is 
determined, the general inverse process of going from toric data to gauge theory data, 



is highly non-unique and a classification of all possible theories having the same toric 
description would be interesting^ Indeed, by the very structure of our algorithm, in 
immediately appealing to the partial resolution of gauge theories on 7L n x 7L n orbifolds 
which are well-studied, we have granted ourselves enough extraneous information to 
determine a unique Qt and hence the ability to proceed with ease (this was the very 
reason for our devising the algorithm). 

However, generically we do not have any such luxury. At the end of subsection 
3.1, we have already mentioned two types of ambiguities in the inverse problem. Let 
us refresh our minds. They were (A) the F-D ambiguity which is the inability to 
decide, simply by observing the toric diagram, which rows of the charge matrix Q t 
are D-terms and which are F-terms and (B) the repetition ambiguity which is 
the inability to decide which columns of G t to repeat once having read the vectors 
from the toric diagram. Other ambiguities exist, such as in each time when we 
compute nullspaces, but we shall here discuss to how ambiguities (A) and (B) manifest 
themselves and provide examples of vastly different gauge theories having the same 
toric description. There is another point which we wish to emphasise: as mentioned 
at the end of subsection 3.1, the resolution method guarantees, upon careful tuning 
of the Fl-parametres, that the resulting gauge theory does originate from the world- 
volume of a D-brane probe. Now of course, by taking liberties with experimentation 
of these ambiguities we are no longer protected by physicality and in general the 
theories no longer live on the D-brane. It would be a truly interesting exercise to 
check which of these different theories do. 

F-D Ambiguity: First, we demonstrate type (A) by returning to our old friend the 
SPP whose charge matrix we had earlier presented. Now we write the same matrix 
without specifying the Fl-parametres: 

/ 1 -11 o -1 o \ 
Q t = -l l o i o-i 

\-l -1 1 1 / 

We could apply the last steps of our algorithm to this matrix as follows. 

(a) If we treat the first row as Q (the F-terms) and the second and third as V ■ U 
(the D-terms) we obtain the gauge theory as discussed in subsection 3.3 and in 

(b) If we treat the second row as Q and first with the third as V ■ U, we obtain 

/-l 1 -1 1 o \ 

d — ( i o o i -2-i which is an exotic theory indeed with a field Ops) 
\oo-ioiiy 

charged under three gauge groups. 

Let us digress a moment to address the stringency of the requirements upon 
13 We thank R. Plesser for pointing this issue out to us. 



matter contents. By the very nature of finite group representations, any orb- 
ifold theory must give rise to only adjoints and bi-fundamentals because its 
matter content is encodable by an adjacency matrix due to tensors of repre- 
sentations of finite groups. The corresponding incidence matrix d, has (a) only 
and ±1 entries specifying the particular bi-fundamentals and (b) has each 
column containing precisely one 1, one —1 and with the remaining entries 0. 
However more exotic matter contents could arise from more generic toric singu- 
larities, such as fields charged under 3 or more gauge group factors; these would 
then have d matrices with conditions (a) and (b) relaxed^. Such exotic quivers 
(if we could even call them quivers still) would give interesting enrichment to 
those well-classified families as discussed in 



Moreover we must check the anomaly cancellation conditions. These could be 
rather involved; even though for £/(l) theories they are a little simpler, we still 
need to check trace anomalies and cubic anomalies. In a trace-anomaly-free 
theory, for each node in the quiver, the number of incoming arrows must equal 
the number of outgoing (this is true for a £7(1) theory which is what toric va- 
rieties provide; for a discussion on this see e.g. |§). In matrix language this 
means that each row of d must sum to 0. 

Now for a theory with only bi-fundamental matter with ±1 charges, since 
(±1) 3 = ±1, the cubic is equal to the trace anamaly; therefore for these theo- 
ries we need only check the above row-condition for d. For more exotic matter 
content, which we shall meet later, we do need to perform an independent cubic- 
anomaly check. 

Now for the above d, the second row does not sum to zero and whence we 
do unfortunately have a problematic anomalous theory. Let us push on to see 
whether we have better luck in the following. 

(c) Treating row 3 as the F-terms and the other two as the D-terms gives 

/ o -1 1 -1 l o \ 
d — I o i o i -2-i which has the same anomaly problem as the one 
\oo-ioiiy 

above. 

(d) Now let rows 1 and 2 as the F-terms and the 3rd, as the D-terms, we obtain 

(Xi X2 X3 X4 x 5 \ 
1 1 -1 -1 , which is a perfectly reasonable matter content. Inte- 
0-1-111/ 

(1 1 o\ 
10010 gives the superpotential W = <f)(XiX 2 X 5 — X 3 X 4 ) 
00101/ 

for some field of charge (0, 0) (which could be an adjoint for example; note 



14 Note that we still require that each column sums to so as to be able to factor out an overall 
U(l). 
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Theory (a) Theory (d) 

Figure 6: The vastly different matter contents of theories (a) and (d), both anomaly 
free and flow to the toric diagram of the suspended pinched point in the IR. 

however that we can not use X\ even though it has charge (0, 0) for otherwise 
the F-terms would be altered). This theory is perfectly legitimate. We com- 
pare the quiver diagrams of theories (a) (which we recall from Figure |3|) and 
this present example in Figure || As a check, let us define the gauge invariant 
quantities: a = X2X4, b = X2X5, c = X3X4, d = X3X5 and e = X\. Then we 
have the algebraic relations ad = be and eb = c, from which we immediately 
obtain ad = eb 2 , precisely the equation for the SPP. 

(e) As a permutation on the above, treating rows 1 and 3 as the F-terms gives a 
theory equivalent thereto. 



(f) Furthermore, we could let rows 2 and 3 be Q giving us d = 




but this again gives an anomalous matter content. 

(g) Finally, though we cannot treat all rows as F-terms, we can however treat 
all of them as D-terms in which Q t is simply A as remarked at the end of 

Section 2 before the flow chart. In this case we have the matter content 

(1-11 -1 \ 
"j -1 1 "i* I which clearly is both trace-anomaly free (each 
1 0-10 / 
row adds to zero) and cubic-anomaly-free (the cube-sum of the each row is also 

zero). The superpotential, by our very choice, is of course zero. Thus we have a 
perfectly legitimate theory without superpotential but with an exotic field (the 
first column) charged under 4 gauge groups. 



We see therefore, from our list of examples above, that for the simple case of the 
SPP we have 3 rather different theories (a,d,g) with contrasting matter content and 
superpotential which share the same toric description. 



Repetition Ambiguity: As a further illustration, let us give one example of type 
(B) ambiguity. First let us eliminate all repetitive columns from the Gt of SPP, giving 
us: 

f\ -1 l\ 

G t — i i o i o , 
\i i i i xj 

which is perfectly allowed and consistent with Figure |[ Of course many more pos- 
sibilities for repeats are allowed and we could redo the following analyses for each of 
them. As the nullspace of our present choice of G t , we find Q t , and we choose, in 
light of the foregoing discussion, the first row to represent the D-term: 



Qt - ( 1 - 2 



-l o l c 

10 



Thus equipped, we immediately retrieve, using our algorithm, 



1000 0\ / Q 1' 

i , , , ,-, I XT* 10 2 1 1 rp _ I -1 1 

1 1 " 1 K = 1 o i o o o 1 - 1 

00111/ V 2 1 0, 



Xi Xi X3 X4, x 5 
1-11-10 
-11-11 



We see that d passes our anomaly test, with the same bi-fundamental matter content 

as theory (d). The superpotential can be read easily from K (since there is only 

one relation) as W = </>(A| — A 3 X 4 ). As a check, let us define the gauge invariant 

quantities: a = XiX 2 , b = c = X 3 X 2 , d = A 3 A 4 and e = X 5 . These 

have among themselves the algebraic relations ad = be and e 2 = d, from which we 

immediately obtain be = ae 2 , the equation for the SPP. Hence we have yet another 

interesting anomaly free theory, which together with our theories (a), (d) and (g) 

above, shares the toric description of the SPP. 

Finally, let us indulge in one more demonstration. Now let us treat both rows of 

our Q t as D-terms, whereby giving a theory with no superpotential and the exotic 
/ -1 1 -1 1 \ 

matter content d — { 1 -2 1 with a field (column 2) charged under 3 
\ o 1 1 -1 -l J 

gauge groups. Indeed though the rows sum to and trace-anomaly is avoided, the 
cube-sum of the second row gives l 3 + l 3 + (— 2) 3 = —6 and we do have a cubic 
anomaly. 

In summary, we have an interesting phenomenon indeed! Taking so immediate an 
advantage of the ambiguities in the above has already produced quite a few examples 
of vastly different gauge theories flowing in the IR to the same universality class by 
having their moduli spaces identical. The vigilant reader may raise two issues. First, 
as mentioned earlier, one may take the pains to check whether these theories do indeed 
live on a D-brane. Necessary conditions such as that the theories may be obtained 
from an M = 4 theory must be satisfied. Second, the matching of moduli spaces may 
not seem so strong since they are on a classical level. However, since we are dealing 
with product U(l) gauge groups (which is what toric geometry is capable to dealing 



with so far), the classical moduli receive no quantum corrections^. Therefore the 
matching of the moduli for these various theories do persist to the quantum regime, 
which hints at some kind of "duality" in the field theory. We shall call such a duality 
toric duality. It would be interesting to investigate how, with non-Abelian versions 
of the theory (either by brane setups or stacks of D-brane probes), this toric duality 
may be extended. 

6 Conclusions and Prospects 

The study of resolution of toric singularities by D-branes is by now standard. In the 
concatenation of the F-terms and D-terms from the world volume gauge theory of a 
single D-brane at the singularity, the moduli space could be captured by the algebraic 
data of the toric variety. However, unlike the orbifold theories, the inverse problem 
where specifying the structure of the singularity specifies the physical theory has not 
yet been addressed in detail. 

We recognise that in contrast with D-brane probing orbifolds, where knowing the 
group structure and its space-time action uniquely dictates the matter content and 
superpotential, such flexibility is not shared by generic toric varieties due to the highly 
non-unique nature of the inverse problem. It has been the purpose and main content 
of the current writing to device an algorithm which constructs the matter content 
(the incidence matrix d) and the interaction (the F-term matrix K) of a well-behaved 
gauge theory given the toric diagram D of the singularity at hand. 

By embedding D into the Abelian orbifold C /(Z n ) fe_1 and performing the stan- 
dard partial resolution techniques, we have investigated how the induced action upon 
the charge matrices corresponding to the toric data of the latter gives us a convenient 
charge matrix for D and have constructed a programmatic methodology to extract 
the matter content and superpotential of one D-brane world volume gauge theory 
probing D. The theory we construct, having its origin from an orbifold, is nicely 
behaved in that it is anomaly free, with bi-fundamentals only and well-defined super- 
potentials. As illustrations we have tabulated the results for all the toric del Pezzo 
surfaces and the zeroth Hirzebruch surface. 

Directions of further work are immediately clear to us. From the patterns emerging 
from del Pezzo surfaces to 3, we could speculate the physics of higher (non-toric) 
del Pezzo cases. For example, we expect del Pezzo n to have n + 3 gauge groups. 
Moreover, we could attempt to fathom how our resolution techniques translate as 
Higgsing in brane setups, perhaps with recourse to diamonds, and realise the various 
theories on toric varieties as brane configurations. 
15 We thank K. Intriligator for pointing this out. 



Indeed, as mentioned, the inverse problem is highly non-unique; we could presum- 
ably attempt to classify all the different theories sharing the same toric singularity 
as their moduli space. In light of this, we have addressed two types of ambiguity: 
that in having multiple fields assigned to the same node in the toric diagram and 
that of distinguishing the F-terms and D-terms in the charge matrix. In particular 
we have turned this ambiguity to a matter of interest and have shown, using our 
algorithm, how vastly different theories, some with quite exotic matter content, may 
have the same toric description. This commonality would correspond to a duality 
wherein different gauge theories flow to the same universality class in the IR. We call 
this phenomenon toric duality. It would be interesting indeed how this duality may 
manifest itself as motions of branes in the corresponding setups. Without further ado 
however, let us pause here awhile and leave such investigations to forthcoming work. 



Appendix: Finding the Dual Cone 

Let us be given a convex polytope C, with the edges specifying the faces of which 
given by the matrix M whose columns are the vectors corresponding to these edges. 
Our task is to find the dual cone C of C, or more precisely the matrix N such that 

N* • M > for all entries. 



There is a standard algorithm, given in [27]. Let M be n x p, i.e., there are p n 



dimensional vectors spanning C. We note of course that p > n for convexity. Out 
of the p vectors, we choose n — 1. This gives us an n x (n — 1) matrix of co-rank 1, 
whence we can extract a 1-dimensional null-space (as indeed the initial p vectors are 
all linearly independent) described by a single vector u. 

Next we check the dot product of u with the remaining p — (n — 1) vectors. If all 
the dot products are positive we keep u, and if all are negative, we keep —u, otherwise 
we discard it. 

We then select another n — 1 vectors and repeat the above until all combinations 
are exhausted. The set of vectors we have kept, w's or —it's then form the columns 
of N and span the dual cone C. 

We note that this is a very computationally intensive algorithm, the number of 

I V \ 

steps of which depends on which grows exponentially. 

\n-lj 

A subtle point to remark. In light of what we discussed in a footnote in the paper 
on the difference between M + = M R a and M' + , here we have computed the dual 
of a. We must ensure that !Z + -independent lattice points inside the cones be not 
missed. 
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